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Abstract 

We calculate the cross sections of the baryon and lepton number violating pro- 
cesses. The proper valley method enables us to calculate the multi-boson processes, 
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1 Introduction 



It is well known that the baryon and lepton numbers are not exactly conserved in the standard 
model, due to the chiral anomaly The topologically non-trivial configuration, especially, 
instanton induces the baryon and lepton number violating process as the following; 



where W represents gauge boson, W ± , Z and H represents the Higgs particle. We denote 
the number of generations by rif and the numbers of the bosons by n w , n^. As 't Hooft first 
indicated!], this process is highly suppressed by the factor exp(— 167r 2 /g 2 ) at low energies. 
This is because the process is accompanied by quantum tunneling through the potential barrier 
that separates the topologically distinct vacua. The top of the potential, which is a saddle point 
in the functional space, corresponds to the static unstable solution, so-called sphaleron0, and 
the height of the potential barrier is given by the energy of the sphaleron: E sp h ~ m w /a. In 
naive picture, the tunneling rate grows up as the energy of the initial state becomes higher, and 
it enhances when the initial state has high energies comparable to E sp h. Since E sp h is O(10) 
TeV, the suppression is expected to decrease or vanish and the baryon and lepton number 
violating processes be observable in TeV region. 

In this viewpoint, the high energy behavior of the process ( [1 . 1|) has been studied in recent 
years 0. By the sphaleron-based calculation 0], it was suggested that the scattering pro- 
cesses involved with many gauge bosons and Higgs particles are not exponentially suppressed. 
Ringwald|| and Espinosa||] calculated the baryon and lepton number violating scattering am- 
plitudes semi-quantitatively by the constrained instanton method]?]], and they also suggested 
that the amplitudes of the multi-boson process increase with the energy and reach lpb for a 
center of mass energy of O(10) TeV. Unfortunately, the high-energy behavior of the process 
with many bosons, which is most interesting, is out of the validity of their calculations. One of 
the problem is that they used the constrained instanton. The constrained instanton method 
is not available for the large-radius instanton, which plays an important role in multi-boson 
process. It is pointed out in Ref.|| that the approximation breaks for n w + ~ 4 in the 
constrained instanton method. 

Another problem of Ref. [Q is that the orientation effect of the instanton is not taken into 
account. The integration of the orientation in the Lorentz SU(2)/j rotations is not trivial but 
very important especially for processes with many gauge bosons. Moreover, in considering 
the multi-boson processes at high energies, masses of the gauge boson, Higgs particle and top 
quark cannot be neglected. All the effect remains unevaluated. 



1 + <1 —> (3n/ — 2)q + n f l + n w W + n h R, 
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In this paper, we overcome these problems and evaluate the cross section. Instead of 
the constrained instanton method, we use the proper valley method , which enables us to 
calculate the multi-boson process for n w + ~ 40 || . The integrations of the orientation 
and phase space integration are performed analytically in the extremely relativistic limit and 
non-relativistic limit, and performed numerically in the general cases. The mass corrections 
of the gauge boson, Higgs particle and top quark are taken into account. 

In Section 2, we consider the valley instanton, which is constructed by the proper valley 
method. This configuration gives an important contribution to the path integral in topo- 
logically non-trivial sector of the functional space. In Section 3, we introduce the collective 
coordinates for the valley instanton. In Section 4, we calculate the multi-point Green function 
(ip 4n f A™ w H Uh ) by the proper valley method. In performing the group integration and phase 
space integration, we obtain the cross sections of the baryon and lepton number violating 
processes. In Section 5, we give conclusion and discussion. 



2 Valley instanton in the electroweak theory 

We consider the standard model with rif generations of quarks and leptons. For simplicity, 
we set SU(3) C and U(l) y coupling constants to zero, and set the Kobayashi-Maskawa matrix 
to the identity. Thus we treat the simplified standard model that reduces to an SU(2) gauge- 
Higgs system with fermionsH |§. As for quarks, we have the left-handed doublets q l £ and the 
right-handed singlets q^, where a = 1, 2, 3 is the color index and i — 1, • • • , n/ is the index of 
the generations. Similarly, we have the left-handed doublets l l L and the right-handed singlets 
l l R for leptons. The bosonic part of the action is given by 

S 9= A [ d 4 x trF^F^ (2.1) 



Sh A 



2g 2 

1 J d'x [{D,H)^ {D,H) + i (fftff - v 2 ) 2 } , (2.2) 

where F^ u = d fM A u — d u A^ — i [A^, A u ] and _D M = <9 M — iA^. The masses of the gauge boson 
and the Higgs boson are 

If 1 

The fermionic part of the action is 

S g = J ^xiiqla^D^qL - iu^a^d^UR - id^a^dn (2.4) 
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+y u q[eH*u R - y d q ] L Hd R + h.c), 
Si = J d^xiil^a^D^L - ie R a^e R - yil\He R + h.c), (2.5) 

where e alS = — e a p = ia 2 , = (cr,i) and = (a, —i). The masses of fermions are given by 
m u,d,i = Vu,d,iV- The total action is given by 

S = S g + S h + Y / (s l + j2S q )- (2-6) 

i=l \ a=l I 

In order to calculate the Green function, (ip 4nf A™ w H nh ) , we must construct a series of the 
configurations, which give a dominant contribution to the path integral in the topologically 
non-trivial functional subspace. They are not solutions of the equation of motion. The vacuum 
expectation value of the Higgs field breaks the classical scale invariance and therefore, except 
the zero-radius instanton, there exists no solution of the equation of motion in this model. 



The dominating configurations make a valley of the action [|1(], pj| . Therefore, the proper 
valley method is suitable to construct them ||. The construction of this series of configurations 
is carried out in Ref. 0. We dub these configurations valley instanton. The valley instanton 
is the solution of a generalized equation of motion, the proper valley equation which is called 
in Ref. || by the name of the new valley equation. The proper valley equation for this system 
is given by, 



F a + -±4—F H + Jr^F^ = X e F, 



8Ap8A v v 5A^SW 5A^5H 



A' : 



j:2 c x2 q x2q 



SW5A^ M 5W5H 5W5W 

where the integration over the space-time is implicit. The valley is parameterized by the 
eigenvalue A e that is identified with the zero mode corresponding to the scale invariance at 
v -> 0. 

To simplify the equation, we adopt the following ansatz; 

Ap(x) = ^f- ■ 2a(r), H(x) = v (1 - h(r)) V , (2.8) 



where rj is a constant isospinor, and a and h are real dimensionless functions of dimensionless 
variable r, which is defined by r = \fx 2 j p. The matrix a^ v is defined, according to the 
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conventions of Ref.0, as a pu = r] apv cr a /2. We have introduced the scaling parameter p. The 
tensor structure in (|2.8|) is the same as that of the instanton in the singular gauge [1]]. 

Inserting this ansatz to ( [2.7|) , the structure of and F H ^ is determined as the following; 

= " ^/ a W, F H \x) = ~f\r) V . (2.9) 



At pv — > 0, (|2.7| ) reduces to the equation of motion and the equation for the zero- mode 



fluctuation around the instanton solution. The solution of this set of equations is the following; 

1 / 2 \ V 2 

1 , I r A \ 



a o — ~. ; — o) "<o 



1 + r2 ' " ' VI + ^ 2 / (2.10) 
fa 2Cr* h = Cr 

JO ( 1 + r 2)2' ^0 (l + r 2)3/2' 

where C is an arbitrary function of pv. Note that ao is the instanton solution in the singular 
gauge and h is the Higgs configuration in the instanton background ]T||. We have adjusted the 
scaling parameter p so that the radius of the instanton solution is unity. The mode solutions 
/q and /q are obtained from da /dp and dh /dp, respectively. 

We construct the valley instanton analytically for pv<l. When pv = 0, it is given by 
the ordinary instanton configuration ao, ho, /q an d fo- When pv is small but not zero, it 
is expected that small pv corrections appear in the solution. On the other hand, at large 
distance from the core of the valley instanton, this solution is expected to decay exponentially, 
because gauge boson and Higgs boson are massive. Therefore, the solution is similar to the 
instanton near the origin and decays exponentially in the asymptotic region. We solve the 
valley equation in both regions and analyze the connection in the intermediate region. In this 
manner we find the solution. 

In the asymptotic region, a, h, f a and f h become small and the valley equation can be 
linearized. The solution of equations is 

a{r) = C ir ^-G pmw (r) + -r(r), 
ar v 

h(r) = C 2 G pmh (r) + ±-f h (r), 



(2.11) 

f h (r)=C 4 G mi (r), 
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where Cj are arbitrary functions of pv, v is defined as A e = v 2 vj\ and p Wth are defined as 
Pw,h — m w,hV^- — 2z/. The function G m {r) is 

CCD - (2-12) 

where i^i is a modified Bessel function. As was expected above, these solutions decay expo- 
nentially at infinity. By matching two solutions ( |2.10| ) and ( |2.11 ) in the intermediated region, 
we obtain C\ + C 3 /v = -In 2 , C 2 + C 4 /v = 2n 2 , C 3 = -4tt 2 C and C 4 = 4n 2 C. By more detail 
analysis, we find C,u — 1/4 at pv — > and we observe C, v ~ 1/4 for /w < 1 numerically [§]. 
Inserting the valley instanton to the action, we obtain, 

S= 8 -4+ 2 -^(pv) 2 + ---. (2.13) 
The fermionic zero mode around the valley instanton is given by, 



~x v af (-e rs 77 ts 7^ L (r) + Vr e^d L {r)) , (2.14) 



QlA x > = - 

u Ra {x) = ^-—u R (r)r] a , (2.15) 

111 p 

d Ra {x) = ^-—d R {r)e (2.16) 

Z7T p 

where the Greek and the Roman letters denote indices of spinor and isospinor, respectively. 
For pv < 1, the solution is given approximately, by the following; 

»Lir)={ ^ \ l (2-17) 

-2tt 2 -— G pmu (r), if r>(pm u )" 1 / 2 ; 
r ar 

— : if r <C (pm u ) -1 / 2 ■ 



«*(»•) = < ^ + 1' ' (2.18) 

2n 2 G pmu (r), if r>(pm u )" 1 / 2 , 

and c?l and d R are obtained by replacing m u with m^. The above approximative behavior is 
correct even for pv ~ 1, because the valley instanton keeps its shape unchanged for pv ~ 1. 
We obtain the zero mode for leptonic sector in the similar manner. 



3 Collective coordinates of the valley instanton 

In this section, we give the collective coordinate method for the valley instanton. Because the 
valley instanton breaks the invariance under the translations and the Lorentz SU(2)r rotations, 
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we must introduce the collective coordinates to restore itQ. However, as the valley instanton 
is not a classical solution except the zero radius one, the zero modes for the translations and 
the Lorentz SU(2)r rotations do not exist in general. Therefore the conventional collective 
coordinate method can not be used. The basic strategy of our collective coordinate method 
is to treat the translations and the Lorentz SU(2)r as a kind of "gauge symmetry" as well 
as the original gauge symmetry |13|]. This guarantees automatically the manifest invariance 
under the translations and the Lorentz rotations through the following deformation and do 
not require the zero modes. 

In addition to these collective coordinates, the valley direction, namely the p integral, must 
be brought into. This is because along the valley direction the action of the system varies 
gently. The corresponding integration can not be approximated by the Gaussian integration 
and must be performed more exactly. The total number of the collective coordinates is eight 
(four collective coordinates from the translations, three from the Lorentz SU(2)r rotations and 
one from the valley direction) and the same as that of the BPST instanton. 

First we rewrite the gauge field and the Higgs fields as 

A^x) = A^(x;p) +Q ll (x;p), (3.1) 
H(x) = H(x;p) + Q H (x;p). (3.2) 

Here A^x; p) and H(x;p) are the valley instanton and Q^(x\p) and Qh{x',p) are quantum 
fluctuations around them. As was shown in (|2.8|) , the tensor structures of the valley instanton 
are 

A,(x;p) = ^p-2a(r), (3.3) 
H(x;p) = v(l - h(r))rj, (3.4) 

where r is r = \/x^/ p. For the fermions, we consider them as the quantum fluctuations. 

Substituting (|3.1| ) and ( |3.2|) for the original Lagrangian, we obtain a Lagrangian for the 
quantum fluctuations. Due to the invariance of the original Lagrangian, it has the invariance 
under the gauge rotations, the translations and the Lorentz SU(2)r rotations. The transfor- 
mation law of the quantum fluctuations are as follows. 



1 In Higgs phase, the collective coordinates for the global SU(2) gauge rotations are not needed as they are 
not symmetries. |l2| 
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1. gauge rotations: g(x) 



u g R ( x ) = u R (x), v%(x) = u R (x), 
d 9 R (x) = d R (x), djf(x) = d R (x), 

Ql( x ) = 9(x)q L (x), q[ 9 (x) = q{(x)g\x), (3.5) 

Qfa; p) = ~ id »9 {x)g\x) + g(x)( A n(x; p) + Qni x \ p))g\ x ) - A »( x ; p), 

Qh( x ', P) = g( x ) ( H (x; p) + Q H (x; p)) - H(x; p). 



2. translations: x 



u x R °(x) = u R (x + x ), Ur°(x) = u R (x + x ), 
d x R °(x) = d R (x + x ), d%°(x) = d R (x + x ), 
Ql ( x ) = q L {x + x ), q[ x ° (x) =q[{x + x ), 
Q x n°( x ', P) = Qv( x + ^o; P) + A^x + x ; p) - A^x; p), 
Qh{x; p) = Q H (x + x ; p) + H(x + x ; p) - H(x; p). 



(3.6) 



3. Lorentz SU(2) R rotations: U 

= U^iR-'x), <(*) = u^R-'x), 
d U Ra {x) = U?d R p(R- l x), dZ(x) = d Ra (R- l x), 

Q U A X ) = q2(R- lx ): = q?(RT 1 x)Ul a , (3.7) 

Q"(x; p) = ^Ua^lf ■ 2a{r) + R^Q^x; p) - A^x; p), 

Q U h( x ;p)=Qh(R' 1 x; P ), 
where R~ l x = R~lx v and the SO (A) matrix R is defined as follows; 

RtivO v = a^U, (3.8) 
and we have used the following equations. 

R^A^R-'x; p) = ^rfa^U ■ 2a(r), H(R' l x; p) = H(x; p). (3.9) 

The transformation law of II and e R is the same as that of qi and d R respectively. 

We will treat these transformations as gauge symmetries. To extract the "gauge orbits" 
from the path integral, we use the Faddeev-Popov technique. For simplicity, we abbreviate the 
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whole original fields by the single filed the valley instanton by $o an d the whole quantum 
fluctuations by Q. We define the Faddeev- Popov determinant as 



A FP (Q)- 1 = J VgdUdx dp5 (F(Q^ u (x;p))) 5 (j d 4 x^- 



5S 



N(p) 5$ (x;p) 



Q**° u (x;p)), (3.10) 



where Q 9X ° u (x; p) is the quantum fluctuations transformed successively under the gauge rota- 
tions g(x), the translations x and the Lorentz SU(2) R rotations U, and 



5 (F(Q^ u (x; p))) = J] S (F„(Q^ u (x; p))) \{ 5 (F^ s (Q^ u (p)) 

x v=\ 

i\s{F^{Q 9xou {p)) 

/ \ 2 

N(p) = I d A x <),S 



0=1 



d x 



1 / 5S 



5S 



5S 



2 \5A lx (x;p) i 



SW(x;p) J \5H(x;p) 



(3.11) 



5S 



Q 9XoU (x ]P ) 



dS QT u ^p) + t^-,Q 9 h oU ^p) 



5S 



5H(x;p) 



Qt oU (x;p)- 



5W(x;p) 

Here F ga _ uge , F Lorentz and -F trans are gauge fixing functions, for example, which are given by 



F gange (Q(x ]P )) = (D;Q^y(x;p), 

Kaaa(Q(j>)) = J d ^ X tr {fp( x \ P)Qnfa p)} , 

Fl oientz (Q(p)) = [ d 4 xtr{f;(x;p)Q,(x;p)} 



(3.12) 



where f"(x;p) and f®(x;p) are appropriate functions for the gauge-fixing. The first delta- 
function in the right-hand side of ( |3.1C| ) gives the gauge fixing and the second one restricts 
the quantum fluctuation to be orthogonal to the valley direction. 

Now we consider the following Green function; 

(V>(xi) • • ■ ip(x4n f ) A^iyi) ■ ■ ■ A„ nw {y n J <j>iz x ) ■ ■ ■ <p{z nh )) 

= M f Dd>^i) ■ ■■i>(x 4nf ) • • • A, nw (y n J 4>{ Zl ) ■ ■ ■ 4>{z nh ) e~ s ^\ (3.13) 



where M is an appropriate normalization factor and ip{x) is the Dirac fermion corresponding 
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to quark or lepton, namely, 

{(ur(x), qLi(x)) for quark in the up sector 

(dn,qL2( x )) for quark in the down sector , , 

(0, for lepton in the up sector 

(eii(x), Il2{x)) for lepton in the down sector 

and 4>(x) is the shifted Higgs field, 4>(x) = H(x) — vrj. Inserting the following identity 

1 = J VgdUdx dpA FP (Q)5 (F(Q^ u (x;p))) 6 (j d'x ^ ^—^ Q^u {x . p ^ j (315) 
we obtain 

(^(xi) • ■■ip{x Anf ) A w {yi) ■ --A^ (y nw ) </>(zi) ■ ■■4>{z nh )) 

= N j VgdUdx dp J V<S>A FP (Q)5(F(Q°»° u (x;p)))5^ d 4 x ^ ^ p) Q 9X ° U (x; p)j 
x^(xi) • ■■it>(x 4n/ ) ^(j/i) • ■■A, nw (y nw ) <f>( Zl ) ■ ■■ct>{z nh ) e~ s ^\ (3.16) 
Substituting the inserted fields in the integral for 

ur^) = ul S u% oU (R(xi - x )), u^(xi) = u^ oU (R(xi - x )), 

<M*0 = Ufd 9 R x f{R{x t - x )), d^{ Xi ) = d^ u {R{x t - x )), 

qt(xi) = gHxJtfo^Wxi - x )), q£( Xi ) = q[ 9X ^(R(x t - x ))U$g( Xi ), 

eRp{ Xi ) = Uf ] e 9 R X s )h T (R( Xi - x )), e R ^(xi) = e^ 9XoU (R(xi - x )), (3.17) 

l$( Xi ) = ^(sOJf ™ - a*)), l[ a (xi) = lT° UP {R{xi - x Q ))U^g{ Xi ), 

A M = r 1 — ^T^gHviW^^giyi) ■ 2a(r yiX0 ) +ig' l (y i )d lx g(y i ) 

{yi — x o) 

+gK yi )R-lQr u {R{ yi - x„); pHVi), 
<f>(zi) = 9*( z i) v ( l ~ H r zixo))v ~vr) + g\zi)Q 9 H oU {R{zi - x ); p), 

where r yiXQ 

= \J{Vi - xo) 2 /p, and using the relations 

= S($ + Q 9XoU ), 
A FP (Q) = A FP (Q gxoU ), (3-18) 
£>$ = VQ = VQ 9XoU , 

all the integration variables which appear in the path integral become Q gx ° u . Then we can 
rename it bake to Q. Furthermore, since the physical observables are invariant under the SU (2) 
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gauge rotations, we can simply set g(x) — 1 in the path-integral instead of the integration of 
g(x). This may breaks the translation and Lorentz invariance of the Green functions which are 
not invariant under the £77(2) gauge rotations, however, the lack of the explicit symmetries 
in such quantities is a gauge artifact. For physical observables, the translation and Lorentz 
invariance are not broken. Finally, we obtain 

(^(xi) • • • i>(x 4nf ) A^( yi ) ■ ■ ■ A^ nw (y nw ) 4>(z x ) ■ ■ ■ 4>{z nh )) 
= N j dUdx dp J VQA FP (Q)6(F(Q(x;p)))s(j d A x ^ ™ Q(x;p)\ 



x^(xi) • • • i/>(x4n f ) 7L M1 (yi) • • • A^ w (y nw ) <p( Zl ) 
where the inserted fields in the integral are given by 
( UWuRpWxi-xo)) 

( UWdwiRixi-Xo)) 

\ <lL2( R ( x i ~ x o)) 

( lURixi-xo)) ) 

( U^eRfsiRixi-Xo)) 
{ l* 2 {R( Xi - x )) 



lj}{Xi) 



i> ] {xi) = < 



N(p) 5$o(x;p) 
4>{z nh )e~ s ^\ 

URa(R(Xi ~ X )) 

q%(R(xi-Xo))U$ 

d ] m {R{xi - xo)) 
q? 2 (R(xi-x ))U$ 


lW{R{xi-x Q ))U% 

^Rai R i X i ~ X o)) 

if 2 (R( Xl - x ))m 



xo) 



Ua^W ■ 2a(r yiX0 ) + Rj / Q v (R(y i - x ); p), 



{Vi - xor 

(j)(zi) = -vh(r ZiX0 )rj + Q H (R(zi - x ); p). 



(3.19) 



4 Calculation of the cross section 



In this section, we calculate the cross section for the process of the type q + q — » (3n/ — 2)q + 
nfl + n w W + n/jH, using the valley instanton. From Section 3, the Fourier transform of the 
Green function is given by, 



G(p,q,k) 



87T 2 /g 2 -8-n w ^-n h /2 



dU J dpe- 27T2 ^ 2 / x p 2n f~ 5 F(pv,pij) 



Xip(p 1 )---A fl (q 1 )---(f)(k 1 )--- } 

where G(p,q,k) is defined by, 

(2tt) 4 5 V + • • • + qi + ■ ■ ■ + h + ■ ■ -)G(p, q, k) 

d 4 xi • • • d A y x ■ ■ ■ d 4 Zl ■ ■ ■ e vi-*i+-+<«i-wi+-+**i-*i- 
x(^(x 1 )---A(y 1 )---(t)(z 1 )---). 



(4.1) 



(4.2) 
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The dimensionless function F(pv,pp) includes the contributions from the fermionic determi- 
nant, bosonic determinant and jacobian factor: 

= % I VQAMQ)S {F{Q(r ' P)))S (I dix W)^h) Q(x ' p) ) e ' slt " +m (4 - 3) 

We introduce the renormalization point /i, and g and A in ( |4.1| ) are renormalized coupling 
constants. In the leading order, and <ft in ( |4.1| ) are given by the valley instanton and ip is 
given by the zero mode around the valley instanton. The Fourier transforms of the fermionic 
zero mode, ip(p) is given by, 

2,i P ( -5*/^ ) 1 + . . . , (4.4) 



m 



f 



where • ■ ■ includes the regular term in the limit: p 2 + m? — > 0, and \ a is the constant spinor, 

(0, —1) for if) = u, v and (1,0) for if) = d, e. We have used the notation: j> = p^a^. The Fourier 
transforms of the valley instanton, is given by 

A„ = — vr ip f +■■■. (4.5) 
v q z + m A w 

The Fourier transform of the shifted Higgs field is given by, 

* = -bvwh? + - (4 - 6) 

By (PD-(p|), we obtain 

G{p q k) — (— i) nw ■ 2 An i'~ nh ii An f +2nw+2nh g~ s ~ nw e~ s ^ 2 l g2 v~ &ns ~ 2nw ~ nh+A 

xD(^;X) f dUL(U;p,q,k) (4.7) 
v J 

where D(^; A) and L(U;p, q, k) are defined by, 

D(^; A) = \- nh/2 f rf(pt»)e" 2 " 2( ^ )2/A (pw) 6n / +2n - +2n ' 1 -V- n -- nh F(pt;,p/i), (4.8) 



m P ,q,k) u[ ) P 2 +m 2 n ;| + ^ n fc?+m 2 

We ignore the momentum dependence that goes to zero on mass-shell. 



(4.9) 



Performing the LSZ procedure, we obtain the invariant amplitude T n nw>nh . Summing up 
the polarization and charge state of the gauge field, isospinor and spinor of the fermion field, 
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we obtain the amplitude; 

E\rp |2 _ C2&nj+n w -2n h 8nf+4n w +4n h ( 2\~n w -S -167r 2 /g 2 -2(6n f +2n w +n h ~-4) 



(4.10) 



i=i 



The cross section is given by the following; 



(7 



1 I ft 



n 



n 



njn h ! 2 4 7 ^ (27r) 3 2£ Pi ^ (2^2E qj ^ (27r)32£ fc; 

1 



(27T)V 4 )(p in -p, 



out ) 



X 



4\/(Pi ' Pa) 2 - mfm 2 



1 



2-4n™-67i h -5^2 7r ^- 4ra /+ n «'+ n h+ 1 0^2^-n tl ,-8 e -167r 2 /g 2 w -2(6n/+2n„+n h -4) 



xD 2 (-;A)/ 2n , n „ )% (s) -s \ 



(4.11) 



where s = pf n = (pi + P2) 2 and the masses of fermions of the initial state are ignored. We 
denote the sum of the momentums of the final state by p ut- All the information about phase 
space and group integration is encoded in the function l2n f ,n w ,n h (s). We define h im ,n{ s ) by, 



J ,_ 3 2±L p . j=1 Zh, q . l=1 



Pout) 



2/ 



2 m / dUHtrfaU) II(-)«trM^) 



i=i 



3=1 



n#n§ i ii^ (4, (p.n-pou,) 

j=3 ZJZ/ Pi i=l ^Qj 1=1 

„ 21 m r , 

J „■ 1 ,■ 1 L 



i=l 



3=1 



(4.12) 



To calculate the cross section, we must evaluate two functions, D(/i/v;X) and Ii )m)Tl (s), 
which are given by ( |4.8|) and ( |4.12|) , respectively. The former includes p integral and the latter 
includes the phase space and group integral. 



4.1 Phase space and group integral 



We consider the Laplace transformation of Ii jm ,n{ s ) 



P2 are independent variables, and at last we input p 



14]| . At first, we assume that p- m , pi and 
= Pi + P2- The Laplace transform 
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®i,m,n(ci',Pi,P2) is given by, 

^,m,n(a;Pi,P2) = / d A p m e _a ' Pin /i )m>n (p in ;pi,p 2 ) 



where 



(f) f {a,U) 



d 3 p 
2E P 
d 3 q 



m. 



tr(0U)e- a - p = 27r—t-K 2 (am f )tr(fiU), 



a 



(4.13) 



(4.14) 



m„ 



-a-q 



2nm 3 



cr 



'-K 3 (am w ) 



{tr$U)f 



a 



d 3 k _ a . k 2rrm h 

e = Ki(am h ). 



(4.15) 
(4.16) 



2E k a 

In the above, K v is the modified Bessel function and a = \fo?. Then we can perform the 
group integration easily, 

2 

— a 

2l-2+2i 



J dUtityJJfafaU) [ti{(/xU)) 21 2 {tr(^C/) 

m „ _ _ _ 

= £ m C< / dUtrfaUfafaU) {tr$U)y ^ (-a 2 )™"' 
= {CoiPi ■ a)(p2 • a) + Ci(pi • p 2 )a 2 } a 2 ^ 2 \ 



where 



Co = 4E(-r"-a 



i=0 



Ci = 2E(") 



i.Ci 



(i + i + iy.(i + i - 2y: 

{2(Z + <-!)}! 



(4.17) 

(4.18) 
(4.19) 



Finally, using the inverse Laplace transformation, we obtain the function I^ m , n - We can 
evaluate Ii >m>n analytically in the case of the extremely relativistic limit and the non-relativistic 
limit. In the extremely relativistic limit, l\, m , n is given by, 



J _ /nrER 3Z+2m+n-4 



(4.20) 



where 



^yER _ nr. 

l,m,n 



i) 



2l+m+n-3 



{(31 + 2m + n-3)\(3l + 2m + n- 5)!} _1 



i=0 



m—i /~t 



(UTTWT^T)! f 1 + (' + ! »( 3 ' + 2m + " " 4 » -< 4 - 21 > 
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In the non-relativistic limit, we obtain 

W = C^mfmWrrCfWi - Mf^ m+ ^~ w • (4.22) 

where 

Cf^ n = | (6Z + + - 10) !p • 2 ^+( m +"V 2 - V' +3 < m+ ")/ 2 - 3 

We assume the fermions are massless except for top quark and define M = mf+n w m w + nhmh- 
In the general cases, we can perform the inverse Laplace transformation numerically, using 
the steepest descent approximation. 

4.2 p integral 

We evaluate the function D(p/v;X) defined by (f4.8|) . In (|4.8Q , F and v depends on pv, 
non-trivially. However, according to Ref.||, the valley instanton is nearly equal to BPST 
instanton in the region pv < 1 and therefore we approximate the forms of F and v by that of 
BPST instanton, if p integral is dominated by the region pv < 1. 

For pv < 1, 

v ~ i F(pt;,p/i) ~ c( PyU )( 43 ~ 8 ™/)/ 6 , (4.24) 
where c is a numerical constant |J 

c ~ 861.94e- ' 997n ^2 10+6n ^7r 6+4n /. (4.25) 
We can calculate p integral by ( |4.24 ) and obtain 



// f \ \ / n\ (43-8n f )/6 

D (E. X) = 2 2 «/+2n,-i cA -n h/2 _M ( ^ ) r(t), (4.26) 



2n 2 ) {vj 

where t = 7n//3 + n w + + 19/12. The integral is dominated by the contribution around 
the saddle point; 

^ = 2>} ' < 4 ' 27 ' 

Therefore, we can check the consistency of the approximation by this value of the saddle point. 
The saddle point pv s depends on the number of bosons monotonously. The valley instanton 
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Figure 1: The energy dependence of the cross section for n w = 75 and nh = 0. The solid 
line denotes the result by the saddle point approximation. The dashed line and the dotted 
line denote the approximate result in the non-relativistic limit and extremely relativistic limit, 
respectively We set n f = 3, fi = m w = 80GeV, g 2 = 0.42, m f = 180GeV and m h = 200GeV 
and neglect the other fermion masses. 

is quite similar to the original instanton even at pv ~ 1, and therefore the approximation 
is valid for n w + rih ~ 40, where we assume rrih = m w and use g 2 ~ 0.42. On the other 
hand, the constrained instanton is deviated from the original instanton at pv = 0.5 |§ and so 
the approximation breaks for n w + ~ 4. Therefore we can calculate the amplitude of the 
multi-boson process, using the valley instanton. Finally, we obtain 

r 2 

a = o- w fc~ 7 (27r)~ 40wf/3 ~ 3n "'~ 3n/l+11/3 r(t) 2 
n w \n h \ 

X (g 2 ) ^n^n^ fjL\ (43 -^ )/3 (4 2g) 
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Figure 2: The cross section for the various n w and n h at s 1 / 2 = 15.55 TeV. 
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n w =40 n w =75 /?>, =100 



Figure 3: The energy dependence of the cross section for n w =0, 20, 40, 75 and 100, and Uh—0. 
A dashed line denotes the unitarity bound. 
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4.3 Results 



We show the result in Fig.[l] for the case where n w = 75 and rih = 0. In Fig.|], the solid line 
denotes the numerical result by the steepest descent method. The dashed line and the dotted 
line denote the approximate result in the non-relativistic limit and extremely relativistic limit, 
respectively. We show cross sections for the various values of n w and in Fig.|2| at ^=15.55 
TeV. From this, we understand that the gauge boson plays an important role rather than 
the Higgs particle. In Fig.^, we show the energy dependence of the cross section for n w =0, 
20, 40, 75 and 100, and rih=0. We show the unitarity bound as a dashed line in Fig.[3]. 
The approximation breaks in the region of the energy where the cross section overcomes the 
unitarity bound. In the region, an interaction between the gauge bosons in the final state [§ 
and interaction between multi-instantons [15] play important roles, which are not evaluated 
in our analysis. 



5 Conclusion 

In this paper, we have applied the valley instanton to the baryon and lepton number violating 
processes. The constrained instanton method was used in the past works, but the processes 
with many bosons, which are expected to be observed, are out of the validity of the method. By 
using the proper valley method, we have given a reliable calculation of the process with many 
bosons. First, by treating the Lorentz rotations and the translations as the gauge symmetries, 
we have introduced the collective coordinates for the valley instanton which is not a classical 
solution. Then we have calculated the cross section of the baryon and lepton number violating 
processes. Furthermore, we have taken into account the instanton orientation effect, boson and 
top quark mass effect, which are not considered yet. Then we have obtained more accurate 
result than earlier works. The final result matches the qualitative estimate |§. The cross 
section grows up, as total energy becomes higher, and the processes with many gauge bosons 
and no Higgs particle grow up more rapidly. 

For high-energy processes associated with many bosons, we see the unitarity is broken, 
which means that our analysis is not appropriate in such cases. It is because we don't take 



into account the multi- instanton effects |T5j or the interaction among the gauge bosons. By the 
proper valley method, new phenomena might be revealed in the subjects. One of the reason 
is that there appear new terms in the interaction between the valley instanton and the anti 
valley instanton. Our present study is the first step toward these subjects. 
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Appendix 

A Faddeev-Popov determinant 



If we choose the gauge conditions as (|3.12| ) and ignore higher loop corrections, the Faddeev- 
Popov determinant becomes 



A FP (p) = 8tt 2 • Det 



8tt 2 • Det 



A ab {x,y;p) B av (x;p) C ah {x;p) 

D^ b (y;p) E^(p) F»\p) X»(p) 

G ab (y;p) H au \p) I a \p) Y a {p) 

W{p) 

A ab (x,y;p) B au (x;p) C ab (x;p) 

D^ b (y-p) E^(p) F^ b {p) x \W{p)\, 

G ab {y-p) H av (p) I ab {p) 



where 



A ab (x,y;p) = (D x D° x )j(x-y) 



{d 2 x 6 ab + 2e acb A c x (x; p)d x + e acd e deb A c x (x; p)A e x (x; p)} 8(x - y), 



B av 


(x; p) = 


— 2tr | 


r a a 
% — 
2 


A x (x 


p) 


C ab 


>; p) = 


-2tr | 


" a a 
i — 
. 2 


A\(x\ 


P) 



d v A x (x;p)j , 



D»\y; p) = -tr D°J^y; p) 



a 



E^{p) = / d 4 xti{r x {x-p)d v A x {x-p)}, 



F» b (p)= d'xtAf^x-p) 



.<y 



i—,A x (x;p) 



G ab (y;p) = -tT\Dlr x (y;p) 



a 



(AA) 



(A.2) 
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H au {p) 

r\ P ) - 
x"(p) - 

Y\p) - 
W(p) = 



(fx tr {ft{x; p)d v A x (x; p)} 
d 4 xtr{f a x (x;p) " 



d xtrU^(x;p)- 



d 4 xtr \f a x {x-p) 



dp 
dA x (x;p) 
dp 



N(p) 



d x 



5S dA a x (x;p) SS dH(x;p) + dH^(x;p) 5S 



5A x (x;p) dp SH(x;p) dp 



dp 5W(x;p) 



When pv — > 0, the instanton measure given by our collective coordinate method is the 
same as one given by 't Hooft[|l|]. To see this, it is convenient to take f£(x; p) and f^(x; p) as 



r = d v AJx;p)-D°A v (x;p), 



.a 



-D G 



2 x 2 + p 2 



(A.3) 
(A.4) 



where D° is covariant derivative under the background of the valley instanton; D Q * = d^ * 



— i[Afj,(x; p), *]. Then it holds 

E^{p) = J d i xtr{mx;p)r x (x;p)} + J d*xtr {/j»(x; p)D° x A u (x; p)} 
F»{p) = f d i xti{r x {x-p)f b x {x-p)}+ f d i xti\r x {x-p)Dl 



a b p 2 



2 x 2 + p 2 

H a »(p) = J d i xtr{r x (x;p)r x (x;p)} + J d*x tr {f a x (x; p)D° x A v (x; p)} , 
I a \p) = j d*xtr{r x (x;p)f b x (x;p)}+ j d'xtr \f a x {x; p)D 



o b p 2 



2 x 2 + p 2 



(A.5) 
(A.6) 
(A.7) 
(A.8) 



When pv — > 0, the valley instanton becomes the BPST instanton and f£ and /" become the 
zero modes corresponding to the global gauge rotations and the translations, which satisfy 



Therefore, as pv — > 0, we obtain 

D» b (y; P )^0, G ab (y;p)^0, 



(A.9) 



E^(p) -> / d*xtr{f x (x;p)r x (x;p)} = S^ftf ■ ~ 



1 
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F^(p) - J d 4 xtr{f x (x;p)f b x (x;p)} = 0, 
H °"(p) -> | d 4 a;tr{/ A a (a;;p)/j:(x;p)} = 0, 
r 6 (p) - / d 4 *tr{/ A >;p)/ A W)} = M/aI| 2 • ^. 



(A.10) 



From O and ( pOOl) , it holds 
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<L4 M (x;p) 
as pt> — >■ and 



p-u 2 fc^fap) 
4A <9p 
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6H(x; p) 



pv 2 dH(x; p) 
4A dp ' 



(A.ll) 



d 4 x 



dH(x; p) 



dp 



0{pv) / d A x 



v dp # 



(A.12) 



Then 



\ 



d 4 x 



'^(x;p)' 



<9p 



|W(p)| 

Therefore, we obtain 

Afp(p) -> 8tt 2 • 2~ 7 ■ \DetA ab (x,y;p) 
as pt> — > 0. When we integrate the quantum fluctuations, the following factor appears; 



)■- 


dA» 




dp 


dp 


■ ii/; 



(A.13) 



(A.14) 



.11/, 



Then the Jacobian becomes 



8tt 2 



1 



)' 


dA^ 




dp 



ii/;ii 4 -ii/;ii 3 - 



(A.15) 



(A.16) 



This coincide with one given by 't Hooft, since 

dA 



dp 



4tt, \\r\\^2V2ir, \\f a J^2np, 



(A.17) 



when pv — > 0. Because the determinant parts coming from the path integral agree with that 
of BPST instanton when pv = 0, the instanton measure of the valley instanton coincides that 
of BPST instanton in this limit. 
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